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ABSTRACT 

We study Galactic emission foreground contamination of the OVRO and UCSB 
South Pole 1994 cosmic microwave background anisotropy data by cross-correlating with 
templates of infrared dust emission and new high resolution VTSS and SHASSA 
data. Ha data provide rough upper limits on the level of free-free emission in the data 
sets. The cross-correlation analysis does not contradict a two component foreground 
emission hypothesis, with the two dust-correlated components being free-free emission 
and spinning dust emission. 

Subject headings: cosmology: observation — cosmic microwave background — diffuse 
radiation — dust, extinction 



1. Introduction 

Cosmic microwave background (CMB) observations have begun to provide interesting con- 
straints on cosmological parameters (see, e.g., Ratra et al. 1999a; Rocha et al. 1999; Knox & Page 
2000; Wang & Mathews 2002; Podariu et al. 2001; Netterfield et al. 2002; Pryke et al. 2002; Stom- 
por et al. 2001; Wang, Tegmark, & Zaldariaga 2002; Efstathiou et al. 2002). To determine high 
precision constrains from CMB anisotropy data, it must be carefully cleaned of any contaminating 
foregrounds to extract purely cosmological fluctuations to high accuracy. 
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So far the main sources of diffuse Galactic microwave emission known to be present are ther- 
mal dust, synchrotron, and frcc-frce. The identification and removal of free-free emission is not 
straightforward as it is likely to dominate the other two sources only in a small range of frequen- 
cies (around 60 GHz). Galactic Ha line emission, which is produced by the same warm ionized 
medium (WIM) responsible for free-free emission, can be a good tracer of free-free emission at high 
Galactic latitudes where there is expected to be minimal extinction from dust. Recently it has been 
proposed that there is at least indirect evidence for the presence of another component of Galactic 
emission, spinning dust emission from rotationally excited small dust grains. The model of Draine 
&: Lazarian ( 1998a, b) predicts a microwave emission spectrum for this component that peaks at 
low frequencies (~ 10—30 GHz), with the details of the spectrum, such as the location of the peak, 
dependent on properties of the dust grains and their environment. Specialized observations directed 
towards picking out this component are being conducted and the first tentative detection has been 
reported in Finkbeiner et al. (2002). 

Indirect evidence for this component comes from microwave-infrared (hereafter IR) correla- 
tions. Estimates of the correlation slope a (defined in § 3) between 100 //m dust emission and 
microwave emission obtained from DMR (Kogut et al. 1996a, b), 19 GHz (de Oliveira-Costa et al. 
1998), Saskatoon (de Oliveira-Costa et al. 1997), OVRO (Lcitch et al. 1997), UCSB South Pole 
1994 (Hamilton k Ganga 2001) and Tenerife (de Oliveira-Costa et al. 1999; Mukherjee et al. 2001, 
de Oliveira-Costa et al. 2002) observations are such that they are roughly consistent with a free-free 
emission spectral index over this wide range of frequencies (10—90 GHz) and angular scales (7° 
to 7'). MAX 4 data show some IR correlation at 100 GHz in one field (Lim et al. 1996; Ganga 
et al. 1998). No IR correlations were found in the Python V (Coble et al. 1999) or QMAP (de 
Oliveira-Costa et al. 2000) data. (If the component that correlates with the IR dust continuum has 
similar scale dependence as IR dust emission itself, which has a multipole space power spectrum 
oc Kogut et al. 1996b, then this relation will be independent of scale^.) However, if this cor- 
related emission were really free-free emission, then a similar correlation should exist between Hq, 
and IR emission^. So far iif^-IR correlations have been found to be marginal. Correlation slopes 
obtained (Reynolds et al. 1995; Kogut 1997; McCullough 1997; de-Oliveira-Costa et al. 2002) arc 
a factor of a few lower than what is expected based on the microwave-IR correlation slopes and 
have significant uncertainties at present. However, if the microwave-IR correlations are only in part 



^Bartlett & Amram (1998) point out that such a power spectrum imphes that the true sky variance remains 

roughly constant towards small scales, but estimates of this variance in restricted patches of sky may be "pulled" far 
from the true variance by important correlations on scales of the patch size. 

''The relation between free-free and Ha emission is Tff ^ 7 jjK ( lo^ k )" '''^'"^'^ (t^)^^ ("r^)' wavelength A, 
where Tg is the brightness temperature for free-free emission, is the intensity of Ha emission, and the power 
law index on the gas temperature changes from 0.55 for temperature T < 2.6 x 10* K to 0.85 for T > 2.6 x lO** K, 
so that at high temperatures there is a greater contrast in the levels of free-free and Ha emission. The assumptions 
that go into deriving this relation include a fraction of singly ionized helium of unity, a helium to hydrogen ratio of 
0.08, an approximate expression for the Gaunt factor, and further the relation is valid only for lines of sight in which 
interstellar extinction is negligible (Reynolds & Haffner 2002). 
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due to free-free emission, then this may be acceptable. This ratio may also vary across the sky due 
to the variation in the Hydrogen ionization fraction an d optical extinction along different lines of 
sight (de Oliveira-Costa et al. 2002). See Smoot (1998), Bartlett & Amram (1998), McCuUough 
et al. (1999), Valls-Gabaud (1998), Kogut (1999), and Draine & Lazarian (1999) for reviews of the 
relation between and free-free emission and for discussions of how spinning dust emission might 
be responsible for the "anomalous" IR-correlated emission. 

While the cause of microwave-IR correlations is still under discussion, we present a correlation 
analysis of the OVRO and SP94 microwave data with templates of dust emission including the 
12 /im template, and new Ha templates, to see whether we are able to separate free-free emission 
and detect spinning dust emission. These experiments are sensitive to very different angular scales 
and span a wide frequency range ('--' 15 GHz to 40 GHz) where a considerable contribution from 
spinning dust emission may be expected. The foreground templates are introduced in § 2, the 
method is outlined in § 3, the OVRO foregrounds analysis is presented in § 4, the SP94 foregrounds 
analysis is presented in § 5, and § 6 consists of an overall discussion and conclusions. 

2. The Foreground Templates 

We use the 100 /xm template (Schlegel, Finkbeiner, & Davis 1998) as a tracer of emission from 
large dust grains which are in equilibrium with the interstellar medium and radiate thermally at 
temperatures around 20 K. The 60 and 12 ^um templates (D. Finkbeiner, private communication 
2000) are used as tracers of emission from smaller dust grains that are transiently heated to higher 
temperatures and emit (thermally at these high frequencies and rotationally at microwave frequen- 
cies) to get back to their ground state. That the 12 /xm template is a clean tracer of emission from 
small dust grains and that such emission has a detectable counterpart at microwave frequencies is 
still being tested. Even if the 12 fim template is a clean tracer of the distribution of small dust 
grains, it may not be a good tracer of spinning dust emission at microwave frequencies which de- 
pends not only on the distribution of grains but on various other local factors. We proceed with 
these caveats in mind. It is thought that the distributions of small and large dust grains may be 
similar on large scales, but in some regions of the sky their small scale distribution could possibly 
differ (Weingartner and Draine 2000). Thus, as noted by de Oliveira-Costa et al. (2002), shorter 
wavelength dust template maps might be expected to trace spinning dust at least as well as the 100 
fim map and possibly better on small scales. The templates are labelled ilOO, i60 and il2 in the 
tables that follow. These templates are such that between 100 fim and 60 fim emission decreases 
over 90% of the sky and between 60 /nm and 12 //m the emission increases over 60% of the sky. 
Over all, between 100 fim and 12 fim the emission increases over about 10% of the sky at high 
Galactic latitudes. Emission from large dust grains is expected to drop to insignificant levels at 
12 fim. In spatial structure the 60 nm template is very similar to the 100 fj,m. template over most 
parts of the sky. 

When cross-correlating the OVRO data with dust templates, cross-correlating with the dust 
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column density template (a template of the amount of thermal dust emission expected if the dust 
were all at a reference temperature of 18 K, Schlegel et al. 1998, and is labelled cdilOO in the tables) 
is necessary in order to make the estimates of the correlation slope a (defined in § 3) comparable to 
those from the DMR data say, as the temperature in the NCP region, and especially in the OVRO 
fields, is lower than the global average, resulting in lower dust emission and a correspondingly higher 
inferred correlation coefficient (the factor is 2.7 for OVRO, this was emphasized by Finkbeiner & 
Schlegel 1999). Table 1 lists the dust temperatures and intensities in the OVRO fields, the whole 
of the NCP region (DEC > 81°), and the high Galactic latitude sky {\b\ > 30°). The estimate 
from the Saskatoon analysis may need to be revised similarly. Note also that there is greater dust 
emission in the OVRO fields and this can play an important role in heating diffuse ionized gas. 

The Ha observations for the OVRO region were carried out on October 4, 2000 with the Vir- 
ginia Tech Spectral Line Imaging Camera (SLIC) (Simonetti, Dennison, & Topasna 1996; Dennison, 
Topasna, &; Simonetti 1997). The SLIC uses a 1.7 nm bandpass filter centered on Ha and a dual 
bandpass continuum filter with narrow passbands on either side of the H^ line. Both filters reject 
the [Nil] lines at 654.8 and 658.3 nm. An //1. 2 Noct-Nikkor lens images a 5-degree radius circular 
field on a cyrogenically cooled TK 512 x 512 CCD camera. A separate auto-tracking CCD camera 
provides arcsecond precision tracking. The facility is located at a dark site at the Virginia Tech 
Horton Research Center. 

We obtained seven 6-minute H^ images and six 3-minute continuum images approximately 
centered on the North Celestial Pole (NCP). We also recorded ten continuum fiat field images 
using a fiat field light box described in Simonetti et al. (1996). The images were combined and fiat 
corrected in IRAF^ using standard techniques (Dennison et al. 1997). For off-axis rays the bandpass 
of the interference filter shifts to longer wavelengths limiting the usable field to a 5-degree radius. 
The sensitivity to Ha is quite flat over this field, but was radially corrected using Ha observations 
of a flat Ha light box. Ha calibration was achieved through repeated observations of the bright 
planetary nebula M76, which we have previously determined to have an Ha flux of 7.95 x 10~^^ 
ergs s~^ cm~^. The H^ image was then continuum-corrected through a weighted subtraction of 
the continuum image. This first required offsetting the continuum image by 1.93 pixels to achieve 
registration. (The different filters exhibit slightly different transverse gradients of optical thickness 
and thus displace the images by different amounts.) The Ha image was then smoothed by 0.5 
pixel to match the resolution of the offset continuum image. The continuum-corrected image has 
significantly reduced stellar contamination, although the brightest stars are not well subtracted due 
to saturation effects. 

The Ha line lies within a band of atmospheric OH emission, and because of the bandpass shift 
with incidence angle, our image contains contamination from the OH(6-l)P2(4) line at 656.9nm in 
the center of the image and the OH(6-l)Pi(3) line at 655.4 nm (Chamberlain 1961) in the outer 



®IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science Foundation. 



-5- 



part. The latter emission is manifested as a faint (si 3 R), diffuse ring with peak radius a; 4.°6. An 
empirical model ring derived from azimuthally smoothed observations of a high-Galactic-latitude 
field devoid of bright stars was used to remove the ring in the NCP field. The fit was optimized 
to remove the ring between l.°6 and 3.°0, the region of interest in this study. Some contamination 
remains outside of this range. The bright star, Polaris, may also contribute to the contamination 
in the center of the image. A uniform gradient caused by scattered light as well as Galactic 
starlight was also removed from the image. Because the image contains an unknown, nearly uniform 
background due to geocoronal Ha emission and scattered light, the median value of the final image 
was set equal to zero. 

Coordinates were established using WCSTools (Mink 1997) which was used to fit to stars from 
the Hubble Guide Star Catalog. The final image is calibrated in rayleighs (R), where 1 R = 10®/47r 
photons cm^^ sr~^. 

For the SP94 region the recently released Southern Sky Survey (SHASSA) Ha, observations 
(Gaustad et al. 2001) were used. 



The method of cross-correlating microwave data with available foreground emission templates 
was pioneered by Kogut et al. (1996a,b), Gorski et al. (1996) and Banday et al. (1996) and further 
generalized by de Oliveira-Costa et al. (1997). It assumes that the microwave data consist of a 
super-position of CMB anisotropy, noise, and foreground emission components. 



Here y is a data vector of N pixels, X is an x M element matrix containing M foreground 
emission templates convolved with the experimental beam and processed by the observing strategy 
adopted, and a is a vector of size M that represents the levels at which these foreground templates 
are present in the data (each component of a is thus the correlation slope for the corresponding 
template). The vector n is the noise in the data and xcmb is the CMB anisotropy convolved with 
the beam and processed by the observing strategy of the experiment. The noise and CMB are 
treated as mutually uncorrelated. 

The minimum variance estimate of a, obtained by minimizing = (y — aX)^C~^{y — aX), is 



3. 



The Method 



y = aX + XCMB + n. 



(1) 



a = [X^C-^X]-^X^C-^y, 



(2) 



with errors Aoj 



S ■ • where the matrix S is 



S =< > - < a >2= [X'^C-^X] 



-1 



(3) 



Here C is the total covariance matrix (the sum of the noise covariance matrix and the CMB covari- 
ance matrix). In these equations X and y are actually deviations of the corresponding quantities 
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from the weighted mean with weights given by C~^. The rms amphtude of temperature fluctuations 
in the data that results from the correlation is AT = (a ± Sa)afoj-e, where (Xfore is the rms deviation 
of the corresponding foreground template map. 

Underlying this method of cross-correlating CMB data with maps of Galactic emission are the 
assumptions that the templates perfectly trace the spatial appearance of the respective components 
at the frequencies and angular scales of the CMB data, and that the templates are exhaustive, i.e., 
they are sufficient to explain all the structure in the data. In addition, for deriving the minimum 
variance estimator, we assume that the noise has zero mean and a Gaussian distribution.^ With 
such caveats in mind, this method has been widely used to test for the presence of contaminants 
in CMB anisotropy data. 



4. The OVRO Correlation Analyses 

The OVRO observations consist of 36 fields around the NCP, at a declination of ~ 88°, at 
frequencies of 14.5 and 31.7 GHz, with a double differencing 7' FWHM beam (see Leitch et al. 2000 
for more information on the OVRO experiment). When the observed data are modelled as consisting 
of a CMB component and a single foreground of spectral index ATj(i^) = ATcmb,i + ^rfore,^^^^) 
then given observations at two frequencies, one can solve for the CMB anisotropy component in 
terms of the unknown spectral index of the foreground, 

where a(i^) = x^e^ jie^ — 1)^, with x = hv/kT, is the correction factor for the RJ approximation to 
a black body spectrum. We use this inferred CMB signal and perform a likelihood analysis. Here 
the total covariance matrix is composed of the CMB anisotropy fiat bandpower theory covariance 
matrix, the instrumental noise matrix, and the (one-off diagonal) correlated noise (T„ = 41 jjK. found 
in the 31.7 GHz data (Leitch et. al. 2000). This process assumes Gaussianity of the inferred CMB 
anisotropy signal and the presence of only a single foreground. A constraint matrix to account for 
the removal of a constant offset is inchidcd in the analysis. After marginalizing over (3 and cj„, the 
constraints on the level of CMB anisotropy agree with those found by Leitch et al. (2000). The 
constraints on /3 are not stringent with /3 > — 1.8 ruled out at 2 a. Although the data are unable to 
rule out a very steep spectrum foreground, based on correlating the data with Westerbork Northern 
Sky Survey (WENSS, Rengehnk et al. 1997) 325 MHz maps, Leitch et. al. (1997) conclude that 
the contribution of any (3 < —2.2 spectrum component is negligible. This then leaves a significant 



^As noted in Mukherjee et al. (2001), violation of the assumptions that the templates are perfect and exhaustive 
can be expressed as an additional "missing" template which would by construction be perfectly compensating. If this 
template is uncorrolatcd with the templates used and has zero mean, we would expect only an increase in the noise 
level. If the template is correlated or has non-zero mean the estimate of the correlation slope a would be biased. 
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amount of foreground at 14.5 GHz of ATrms ~ 175 /xK, though it indicates that the 31.7 GHz data 
are almost entirely CMB anisotropy. 

Table 2, in its upper panel, lists results from correlating the OVRO data with each dust 
template individually. Of the three templates of dust emission, the 100 fim template correlates best 
with the 14.5 GHz data, followed closely by the 60 fxm template, and then the 12 fim template, all at 
a high level of significance. The reverse trend may be true for the 31.7 data, however the significance 
of dust correlation with the 31.7 GHz data is small^°. The correlation slope for the cdilOO template 
is a factor of about 3 lower than that for the ilOO template. The results are independent of the input 
CMB power spectrum. The of the fits indicate that the small non-CMB signal in the 31.7 GHz 
data, inferred from eq. (1), is satisfactorily accounted for, and about half of the non-CMB signal 
at 14.5 GHz can be accounted for by the dust templates^^. However, the inferred foreground and 
the fitted foreground have similar spatial structure so it is possible that dust-correlated emission 
might satisfactorily explain all the foreground in the 14.5 GHz OVRO data and the discrepancy 
might be due to uncertainties in the foreground template. If the templates are assumed to be close 
to perfect, then a better job can be done if we are able to cleanly separate the two components of 
emission likely to be present in the data, as discussed below. Another possibility is that there is 
an unaccounted-for foreground, possibly free-free emission uncorrelated with dust. 

The central panel of Table 2 lists results from correlating the OVRO data sets with two 
foreground templates jointly. When correlating jointly with the 12 nm and 100 ^m templates, 
all the correlation in the 14.5 GHz data gets attributed to the 100 nm template. At 31.7 GHz 
there might be slight evidence of the emission getting attributed to the higher frequency template. 
These dust templates have similar spatial structure in the OVRO fields (r = S^-^/(S^^S~.^)°-^, 
the weighted correlation coefficient between the templates indexed i and j, is 0.7, see Figure 1). 
When the 60 and 100 ^um dust templates are correlated jointly with the OVRO data, the high 
correlation between these templates (r = 0.9) results in large uncertainty in separating the two 
dust emission components. (When the 12 jim and 60 /nm templates are jointly correlated the result 
is similar to when the 12 /xm and 100 //m templates are jointly correlated.) 

If the dust-correlated emission were entirely spinning dust emission, then a stronger correlation 
with the higher frequency dust templates should be expected. The results from individual and joint 
template correlations demonstrate that this is not the case for the OVRO 14.5 GHz data. The 12 
)um-correlated emission is subdominant to the 100 /xm-correlated emission in the 14.5 GHz OVRO 



If wc calculate the correlation slope after marginalizing over cr„, the correlation slope so found is well within the 
error bars of the numerical values in Table 2 but with slightly larger error bars. 

^^We have tested for possible biases in this minimum variance analysis method by analyzing 1000 simulations of 
the data. These were made by simulating the CMB anisotropy component from a model and adding noise and a 
certain level of foreground on to it. We find that the weighted mean of the res ulting correlation slopes turns out 
to be what was put in to a high level of accuracy and the weighted rms deviation of the correlation slopes from the 
mean was in good agreement with the error bar for the chosen model. Thus not only is the method not biased, but 
the error bars are also accurate. 
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data, even though the two templates have similar spatial structure. This may happen if the 12 fxm 

template is a better tracer of spinning dust emission, and the 100 fiui template additionally traces 
free-free emission. The numbers in Table 2 would then indicate that there are comparable amounts 
of the two components present in the OVRO data at 14.5 GHz. 

More microwave data should be used to check this deduction, which is only tentative at present 
given the caveats discussed in the earlier section and the shortcomings and possibilities for system- 
atic errors in this correlation method discussed in detail in Mukherjee et al. (2001) and mentioned 
in section § 3. If however this deduction is true then the 14.5 GHz data are consistent with about 
92 ± 16 ijK rms of spinning dust emission and hence, using the 100 /xm-correlated emission results 
for the 14.5 GHz data (Table 2), with roughly 101 it 22 /jK rms of free-free emission at 14.5 GHz. 
This would translate into 19 it 4 rms of free-free emission at 31.7 GHz. Note, however, that the 
uncertainties in estimates of the correlation slope and the corresponding temperature fluctuations 
are comparable at the two frequencies. Hence, since the dust-correlated foreground components 
are present only at low levels at 31.7 GHz the significances of the correlations are low as expected, 
becoming even lower when the data at this frequency are correlated with the 12 /xm and 100 f^m 
templates together. This makes it hard to ascertain which component dominates or even distin- 
guish between two components at 31.7 GHz. However the estimated a values do not contradict 
the above scenario, which would explain why the spectral index of the 100 /nm-correlated emission, 
-2.61° ;q, is more skewed towards steeper values while being consistent with free-free, since the 100 
/xm-correlated emission at 14.5 GHz contains in it a 12 /xm correlated component as well, which 
itself has a spectral index of —1.7~^^'l between the two frequencies, probably consistent with what 
is expected for spinning dust emission. 

When the Ha data are convolved with the OVRO beam and double differenced, the data have 
an rms of 0.83 R in the OVRO fields. To find the visual extinction (which results from a combination 
of scattering and absorption) expected in the OVRO fields, we use the reddening maps provided 
by Schlegel et al. (1998) and find that the average reddening in the 36 OVRO fields is E(B-V) = 
0.297 mag, implying a high visual extinction of 0.921 mag (using R(V) = 3.1).^^ This implies that 
the Ha signal seen is at most io-0-^><0-92i = 0.428 or 43% of the total Ha signal, or that the total 
Ha signal is 2.3 times that observed. This correction for extinction is not accurate as the decrease 
in the Ha surface brightness in a given direction depends not only on the total column density of 
dust but also on the unknown arrangement of the emitting and absorbing material along the line 
of sight (Reynolds & Haffner 2002). Without attempting to make a more accurate correction for 
extinction, 0.83 R would imply an estimated Ha signal of 1.9 R and a free- free signal of 63 /xK rms 
at 14.5 GHz. This level of free-free emission estimated from Ha data is an upper limit because 
some of the Ha signal could be due to incompletely subtracted stars. 



^^The average extinction in the whole of the NOP (DEC > 81°) region is 0.15 mag, in the b > 30° region it is 0.04 

mag, and in the b < —30° region it is 0.05 mag. 

^^If, in an attempt to remove stellar residuals, we set, prior to convolution and double differencing, the signal in 



-9- 



Thus the upper limit on free-free emission inferred from data is less than the free-free like 
100 /im-correlated signal found above, and is much less than the signal in the OVRO data at 14.5 
GHz that was found to have a spectral index like that of free-free emission pointed out at the 
beginning of this section. Could there be a preferential heating of the WIM, given that gas gets 
heated more easily than dust, due to the excess dust emission in this region, such that the 
data could be used to deduce a higher level of free-free emission in this region? For there to be, 
say 100 fiK rms of free-free emission at 14.5 GHz, there needs to be about 3 R of (correlated) 
fluctuations, or only 1.3 R in the presence of 43% visual extinction. If the gas temperature here 
were higher by a factor of 2.3, then the level of if^ fluctuations expected can come down to 0.83 
R.14 

The results from correlating the OVRO data with the Ha template are given in the bottom 
panel of Table 2. When correlating the Ha template with microwave data one may not flnd an 
observable correlation when there is significant visual extinction due to dust grains for which the Ha 
template has not been corrected, and also when the data contains a significant amount of another 
foreground. However, a negative correlation with the H^ template is found when the 14.5 GHz 
data are correlated jointly with the Ha and 100 ^m or 12 yUm templates with somewhat greater 
significance being attributed to the dust correlations than was obtained in the case of individual 
correlations. From this we deduce that the structure in the Ha template partially traces the 
structure in the dust templates, sky-rotations indicate the same. This can also be seen from Figure 
1. Note however that this correlation becomes apparant only after the OVRO double differencing. 
Before the OVRO double-differencing, none of the dust templates show a significant correlation 
with the Ha template. Low //a-dust correlations may be acceptable, however, due to various 
reasons^^, including the possibility that microwave-dust correlations are only in part due to free- 



pixels in the Ha map that show structure greater than 3 a to cr which is the local rms deviation from the local mean, 
where local here is taken to be from half a degree below to half a degree above the OVRO ring, then the rms of the 

resulting Ha signal in the 36 fields drops to 0.39 R (0.48 R before the double differencing). We have repeated the 
analysis with this template and the results and conclusions presented here still hold. 

^''Leitch et al. (1997) conclude that the temperature needs to be higher by a factor ~ 100 in order to explain the 
level of foreground in the OVRO data. This is because they wanted to explain a free-free emission signal of 300 /xK 
(a factor of 3 larger from what is argued for here), they had not taken the high visual extinction in the OVRO fields 
into account (this causes a difference of a factor of 2.3), and had an Ha signal of rms < 0.1 i? (factor of 4 lower 
at least from what is found here), giving in total a difference of a factor of 28 from that found here, which when 
multiplied by 4.2 leads to the factor ~ 100. 

^^lonized gas emits Ha and free-free in direct proportion, with the surface brightness of emission proportional to 
the electron density squared. These emissions only indirectly trace thermal emission from the dust associated with 
the ionized medium which has a surface brightness proportional to the density of dust. Further dust absorbs Ha light 
and so about 25% of high latitude Ha emission is estimated to be scattered light (McCuUough 1997). Lagache et al. 
(2000) have used Ha and HI data to decompose the far infrared thermal dust emission at high Galactic latitude into 
components associated with the WIM and the warm neutral medium (WNM). They confirm that about 20 — 30% of 
thermal dust emission is associated with the WIM. The correlation between ionized gas and dust emission may also 
be angular scale dependent (McCullough 1997). 
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free emission, in patches of the sky with high gas temperature and optical extinction, and on small 
angular scales. Since all these are true of the OVRO fields, estimates of correlations involving Ha 
data in the OVRO fields may be suspect and we do not present numbers for Ha-dnsi correlations 
here. 



5. Foreground Emission in the UCSB South Pole 1994 Data 

The SP94 data were taken in two bands. The Ka band had four channels (centred at 27.25, 
29.75, 32.25, and 34.75 GHz), hereafter K^l, Ka2, Ka3, and i^a4, and the Q band had three (at 
39.15, 41.45, and 43.75 GHz), hereafter Ql, Q2, and QS. The beams had resolution of about 
1° and there are 43 data points in each frcqTicncy channel, all at a declination of —62°, with RA 
ranging from 23° to 67°. The data were obtained using a sinusoidal l.°5 chop with smooth, constant 
declination, constant velocity scans. See Gundersen et al. (1995) and Ganga et al. (1997) for more 
details on the SP94 experiment. Hamilton & Ganga (2001) have reported the detection of 100 
/xm-correlated emission at 1.6 a confidence in the Q band data, resulting in a possible reduction 
of the deduced CMB anisotropy band temperature of Ganga et al. (1997). The Ka band data was 
found to be free of contamination from this template. 

The SP94 cross-correlation analysis here is done ignoring the last 4 SP94 pixels that are close 
to and affected by a set of pixels at RA ~ 69° that have negative values of flux in the 12 ij,m 
map, apparently due to inaccurate source subtraction at that location. Wc show the results of 
correlating the data with the 100 /im and the 12 ^um templates individually because when the data 
are correlated with these templates jointly, then due to the high correlation between the templates 
(r = 0.8, see Figure 2), uncertainties increase further and nothing can be gained, given the already 
low level of contamination in this data. Estimates of the correlation slope are derived by jointly 
correlating all channels of the Ka band, the Q band and also all channels of both bands together, 
following Hamilton & Ganga (2001). This is done because the anisotropy model covariance matrix 
causes non-trivial covariance between pixels of all channels. The noise in a given pixel is correlated 
with noise in the same pixel in other channels of the same band. The noise in the Ka and Q bands 
are uncorrected. We then combine the a values obtained for each template into a representative 
(average) number (a) using both a spectral index n (here a = ao{i'/i'o)'^, where vq is a reference 
frequency) of and the best fit n, as described in Hamilton &: Ganga (2001), for the Ka and Q 
bands separately and for both of them together. The best estimate of a is obtained, in the case 



^^If we include all the SP94 pixels in the analysis, then we get a strong correlation of the SP94 data with the 12 /xm 
template. The pixels in the 12 /jm map that have negative flux values cause a dip in the signal in the last few pixels 
of the SP94 beam convolved template (even when the negative flux values are set to or to a value lower than the 
local mean, assuming that we can not really know the level of diffuse emission here due to the presence of a source 
that has not been accurately removed). There is a similar dip in the microwave signal in the different channels of the 
SP94 data (see Figure 2), and this results in a strong correlation with the 12 /um template. 
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n = 0, from 

_ _ Total(S-ia) 
" " Total(S-i) 

where Total denotes the sum of ah elements of a matrix or vector. This is thus a weighted average 
taking account of correlations between the a values of the different channels. The best fit n is 
obtained by minimizing, over a grid of values of n, = (f* ~ dom)'^T,~^[a — abm), where m is a 
vector with elements rrii = [vi/vg]". 

Results are shown in the Table 3. We see that the Ka band data does not correlate with the 
100 //m template while the Q band data correlates marginally, the best fit spectral index when 
both bands are considered together is positive, and there is a greater than 2 a detection of 100 
/im-correlated emission, in agreement with Hamilton & Ganga (2001). Regarding 12 /im-correlated 
emission, we again find that the Ka band data does not correlate while the Q band data correlates 
marginally, the best fit spectral index when both bands are considered together is positive, and 
there is a greater than 2 a detection of 12 /xm-correlated emission when all channels are considered 
together with the best fit spectral index. The upper limit to the total contamination from dust 
templates in the SP94 data is ATj-ms = 26 it 11 fiK when the results of all channels are combined 
using estimates of the best fit spectral indices. This is an upper limit because the 100 and the 12 
lira correlated emission have been coadded, even though these templates are seen to have similar 
structure (see Figure 2) . We are unable to say from this kind of a correlation analysis which of these 
two components is dominant or even whether these correlations are due to two physically distinct 
components. Most of the contamination in the SP94 data is in its Q band (total dust-correlated 
emission in this band is at most ATj-ms = 32 ± 21 fiK), as a consequence of which the deduced 
SP94 band temperature (Ganga et al. 1997) in the Q band could go down in amplitude by at most 
21%, to 41 it 26 /xK. No clear evidence is found for the contamination of Ka band data. This is in 
agreement with the conclusions of Ganga et al. (1997), also see Ratra et al. (1999b). 

Sky-rotations (correlating the data with template dust maps obtained by rotating the initial 
map around the Galactic poles and by inverting the North and South Poles) indicate that the 
correlations are as significant as the numbers indicate. In other words the correlation that could 
arise from a random alignment between CMB anisotropies and the dust template has been taken 
into account in the analysis by the data and CMB covariance matrices, and the rms of the corre- 
lations from random patches agrees well with the error bars obtained by the correlation method. 
Correlations similar to those obtained using the 100 /xm template are obtained when the data are 
correlated with the column density template for 100 /xm emission, indicating that the dust temper- 
ature in the SP94 pixels is not radically different from the average. The optical extinction in this 
region is also small. 

We analyze the structure of Ha signal in the SP94 pixels using data from SHASSA. These data 
are available at a resolution of 0.°1. No significant correlation is found when the SP94 data are 
correlated with the SP94 beam convolved Ha data. When the SP94 data are correlated jointly with 
the 100 /im and Ha templates, a negative correlation is attributed to the Ha template {Ka + Q: 
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n = 0, a = —3.6 it 4.1 and n = 3, a = —3.9 ± 2.1) while a positive correlation is attributed to 
the 100 template {Ka + Q: n = 0, a = 35.6 ± 34.4 and n = 2.2, a = 72.7 ± 22.8), of greater 
significance than that obtained when the data were correlated with this dust template individually. 
Thus, like in the OVRO case, we find that the structure in the Ha template has something in 
common with structure of dust emission (see Figure 2), and that in such a situation the correlation 
method is unable to separate any iJ^-correlated signal from the rest of the dust-correlated signal. 
The correlation slope between the Ha signal and the 100 /xm dust emission, for example for channel 
Ql, is 0.3 lb 0.2 R(MJy/sr)~^ (10% of the rotated dust patches have greater correlation slope, 8% 
of the dust patches correlated with greater significance). The correlation coefficient is 0.4 it 0.2 (the 
same 8% of the dust patches had greater correlation coefficient than the real dust patch) . Here the 
variance is found by correlating the Ha data with 144 rotated versions of the dust emission template. 
As in the case of OVRO, the Ha and dust data correlate better after the chopping strategy of the 
SP94 experiment has been simulated on these data. More sensible deductions regarding any such 
association can be made by comparing the structure in these maps over larger areas of sky (rather 
than in small Id strips). 



6. Discussion and Conclusion 

Consider the scenario in which dust-correlated emission consists of free-free emission as well 
as another component. The amplitude of dust-correlated emission in DMR data at 53 GHz is 
ATrms = 7.1 lb 1.7 fiK, while the total amount of free-free emission is ATj-ms = 5.2 ± 4.2 /xK using 
multifrequency data (Kogut et al. 1996b). Since at that time no other dust-correlated component 
of emission was expected to be found at these frequencies, it was concluded that the correlated 
component forms at least a third (at 95% confidence) of the total free-free emission and possibly 
even the bulk of it. If instead free- free emission were only a part of the total dust-correlated emission 
(with the 12 /xm and 100 /xm templates being similar in structure at large angular scales) then the 
remaining part is constrained to be 5 ± 4 /xK rms at 53 GHz on DMR scales. If this component is 
independent of angular scale, then this estimate together with the amplitude of 12 |[xm-correlated 
emission at 14.5 GHz (OVRO) and 36 GHz (SP94) implies a spectral index of (3 = -2.21'j';^ between 
14.5 GHz and 53 GHz and /? = -3.0+1? between 36 GHz and 53 GHz (between 14.5 GHz and 31.7 
GHz (3 = — 1.7+5'5 ^ discussed in § 4). These spectral index estimates may be consistent with the 
spinning dust hypothesis, becoming more negative at higher frequencies, although the uncertainties 
in the estimates are significant. For free- free emission, if we use the power spectrum estimated by 
Veeraraghavan &; Davies (1997), Q oc (they estimated the power spectrum directly from 

cleaned and destriped Ha maps and so may have obtained a more reliable estimate of the free-free 
emission power spectrum than did Kogut et al. 1996a), and predict the level of free- free emission 
that the OVRO 14.5 GHz and the SP94 36 GHz experiments should see at their respective angular 
scales given the average estimate of COBE DMR, we get 26 ± 20 /xK rms and 6 it 5 /xK rms for 
OVRO and SP94 respectively. This implies free-free emission of < 66 /txK rms at 95% confidence 
in the 14.5 GHz OVRO data, consistent with the rough upper limits indicated by Ha observations, 



-13- 



though the emission in smaller regions could deviate significantly from this average. We see that it 

is plausible that the 100 template, at least in some regions of the sky, traces free-free emission, 
as well as the structure in the 12 /im template by virtue of the large and small grains being well 
mixed. 

If it is generally true that the 100 /um template correlations are in part free-free, the remaining 
being another component that the 12 fim template traces better (this may hold only in regions of 
the sky where large and small dust grains are well mixed, also our analysis throughout assumes 
that the emission from these grains are well traced by the templates we use which are thus close to 
perfect as discussed), this hypothesis would explain why Veeraraghavan & Davies (1997) obtain a 
considerably lower normalization for the free-free emission power spectrum derived from Ha maps 
than did Kogut et al. (1996a), who used dust-correlated microwave emission (which is only in 
part free-free emission). It would also explain the flatter spectrum of free- free emission that they 
obtained given that the remaining dust-correlated emission makes a contribution to temperature 
fluctuations that seems to be independent of angular scale, from the estimates available to date. 

Figure 3a is a plot of 100 /zm-correlated emission found in microwave data so far. It is similar 
to Fig 4. of Hamilton & Ganga (2001) with the corrected OVRO points with error bars added. 
Spectrally the data do not seem to favour either of the spinning dust or free-free interpretations. 
Note that the set of Tenerife points taken from de Oliveira-Costa et al. (2002), for the region 
|6| > 20° no longer show a rising spectrum, the key signature of spinning dust emission, between 
10 and 15 GHz (though they did find that microwave- i^c correlations were too small to explain the 
microwave-IR correlations in the data). The points taken from Mukherjee et al. (2001) are higher 
because there it was shown that only pixels close to the Galactic plane actually show dust-correlated 
emission and here the estimated a is higher. For such reasons comparisons between experiments 
must be made with caution. 

Figure 36 is a plot of 12 /xm-correlated emission in microwave data. The Tenerife points (taken 
from de Oliveira-Costa et al. 2002) could go up if the estimate is confined to the limited region that 
shows dust-correlated emission at both frequencies (Mukherjee et al. 2001). The figure indicates 
that spinning dust emission might have been detected in microwave data with significant levels of 
emission at about 15 GHz. We however need to check for example whether the correlation slope 
between the 12 /xm template and DMR data lies on the spectral curve shown. In order to correlate 
the DMR data with the 12 /Ltm template we may need to await the cleaning of the template as the 
signal at 12 /im is contaminated significantly by zodiacal emission. Modelling zodiacal emission is 
a complex task and this could take some time. We were able use the 12 ^um template in its present 
state as we consider only small regions of the sky, and in these few pixels far from the ecliptic plane 
there does not seem to be obvious contamination. So, whether the 12 /xm-correlated emission in 
other CMB anisotropy data sets follow the spectral pattern of this figure remains to be determined. 
Higher quality data from the newer CMB experiments, such as CBI and DASI, and that anticipated 
from MAP, should help resolve this issue. 
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Rather than having to use thermal dust emission as an indirect tracer of free-free emission, 
data from new Ha observations are expected to provide a template that more directly traces free- 
free emission. However, converting the Ha data into a reliable template of free-free emission may 
require some care. A careful comparison of structures in a high resolution Ha map, with information 
on extinction and gas temperatures (say from X-ray and ultraviolet observations), to structures 
in microwave and IR data over larger areas of sky is needed to better understand the relevant 
associations. The spinning dust hypothesis itself is being studied by making pointed observations 
towards dust features predicted to emit significant amounts of spinning dust emission (Finkbeiner 
et al. 2002). Such observations can more easily and significantly constrain the spinning dust model. 
Even if these foreground contaminants are present only at low levels in CMB anisotropy data, it is 
important that they be understood and removed to high precision, especially if one wishes to study 
higher order statistics and non-Gaussian signatures in the CMB anisotropy data (see, e.g.. Park 
et al. 1998, 2001; Chiang k Coles 2000; Barreiro & Hobson 2001; Wu et al. 2001; Park, Park, & 
Ratra 2002, Shandarin et al. 2002). 

In summary, we discuss a scenario in which 100 /xm-correlated microwave emission in the 
OVRO and SP94 microwave anisotropy data sets consists of free-free emission as well as another 
component. This is motivated by the facts that free-free emission by itself is not expected to result 
in strong microwave-IR correlations, free-free emission is constrained by Ha emission in both the 
data sets to be less that the amplitude of dust-correlated emission, and that the 100 /xm-correlated 
emission dominates over the 12 /xm-correlated emission in the OVRO 14.5 GHz data. The Ha data 
and dust emission data seem to have some common structure, but dust emission correlates better 
with microwave data. If we use the power spectrum of Veeraraghavan & Davies (1997) for free-free 
emission and assume that the other component is independent of angular scale, as observations to 
date seem to roughly indicate, or if we use the 12 jim template as a tracer of this component, then 
we find that from spectral index considerations this component could be spinning dust emission. 
Other CMB data sets would have to be studied to confirm this. Given that the two components 
may have different power spectra, which component dominates in a given data set would depend 
on both the frequency as well as the relevant angular scales. Higher quality CMB anisotropy data, 
more knowledge about the characteristics of dust emission at high frequencies and at microwave 
frequencies, and high resolution Ha data over large portions of the sky, will all be valuable in 
obtaining a clearer picture of foreground contaminants in low frequency microwave data. 
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Table 1. Dust Temperatures and 100 fim Template Intensities 



Region 


Average Dust Temperature (K) 


il00''(MJy/sr) 


cdil00''(MJy/sr) 


OVRO Fields 


16.2 


5.90 


16.1 


NCP (DEC > 81°) 


16.9 


4.08 


8.41 


|6| > 30° 


17.9 


2.12 


2.56 



^The 100 /Ltm dust emission template. 

'^The 100 fim column density template that assumes a constant dust temperature of 18 K. 



-17- 



Table 2. Results of Correlating the OVRO Data with Foreground Templates 



14.5 GHz 31.7 GHz 



Template a±Sa 


a 

T — 

oa 


AT (uK) 


a ± (5a 


a 

-? — 

oa 


AT (aK) 


cdilOo'^ 


50.8 ± 5.3 


9.6 


134 ± 14.0 


6.6 ±5.3 


1.2 


17.3 ± 14.0 


d 

ilOO 


145 ± 15.7 


9.2 


136 ± 14.8 


19.5 ± 15.2 


1.3 


18.3 ± 14.3 


ieo*" 


1210 ±133 


9.1 


138 ± 15.1 


191 ± 128 


1.5 


21.7 ± 14.6 


f 

il2 


2190 ± 373 


5.9 


91.6 ± 15.6 


578 ± 361 


1.6 


24.1 ± 15.1 


il2 


-246 ± 506 


-0.5 


-10.0 ±20.6 


484 ± 480 


1.0 


20.2 ±20.1 


ilOO 


152 ±21.3 


7.1 


143 ± 20.1 


6.0 ± 20.3 


0.3 


5.6 ±19.1 


i60 


531 ± 336 


1.6 


60.5 ± 38.2 


245 ± 312 


0.8 


27.9 ± 35.5 


ilOO 


87.3 ± 39.6 


2.2 


82.1 ± 37.3 


-7.1 ±37.1 


-0.2 


-6.7 ±34.9 


Ha 


-17.0 ±19.6 


-0.9 


-15.7 ±18.0 


-2.0 ± 17.9 


-0.1 


-1.7 ± 14.8 


Ha 


-36.0 ±19.7 


-1.8 


-29.9 ± 16.4 


-5.3 ±18.0 


-0.3 


-4.4 ± 15.0 


ilOO 


148 ± 15.8 


9.4 


139 ± 14.9 


20.1 ± 15.4 


1.3 


18.9 ± 14.5 


Ha 


-70.5 ±21.1 


-3.3 


-58.7 ±17.6 


-18.6 ±20.0 


-0.9 


-15.5 ±16.6 


il2 


2700 ± 403 


6.7 


110 ± 16.3 


747 ± 404 


1.8 


30.3 ± 16.4 



'^The upper panel in this table refers to individual correlations of the OVRO data with dust 
templates, the central panel considers joint correlations of the data with dust templates and the 
bottom panel refers to correlations of th e data including the Ha template. 

^Correlation slope a in units of /xK(MJy/sr)~^ for microwave-dust correlations, and in units of 
/xK(R)~^ for microwave- if (J correlations. 

^The 100 jim column density template that assumes a constant dust temperature of 18 K. 
'^The 100 /J,m dust emission template. 
^The 60 iJ,m dust emission template. 



The 12 /xm dust emission template. 



Table 3. Results from Combining Channels of the SP94 Data^ 



Ka Q Ka + Q 



Template 


n 


a ± (5a 


a 
5a 


c 

AT 


n 


azkSa 


a 
5a 


c 

AT 


n 


a±Sa 


a 
Sa 


c 

AT 


ilOO 


0.0 


28.9 ± 32.3 


0.9 


8.5 ±11.3 


0.0 


62.6 ± 42.7 


1.5 


23.4 ± 16.1 


0.0 


24.8 ± 34.4 


0.7 


5.0 ± 12.1 




d 








0.8 


63.1 ± 42.3 


1.5 


23.9 ± 16.0 


2.2 


53.4 ± 22.8 


2.3 


20.8 ±8.9 


il2 


0.0 


129 ± 716 


0.2 


-0.2 ± 11.5 


0.0 


1150 ± 799 


1.4 


20.1 ± UA 


0.0 


431 ± 756 


0.6 


0.7 ± 12.2 




d 








3.2 


1190 ± 697 


1.7 


22.1 ± 12.9 


3.0 


854 ± 359 


2.4 


16.2 ±6.8 



^We combine the Ka band channels, the Q band channels, and finally all the channels, using first a spectral index n = and then the best 

fit n for the case when the SP94 data arc correlated with the 100 /im and 12 /./m templates individually. Here Ka band results use a reference 
frequency vq = 30 GHz, Q band results use vq = 42 GHz, and when all channels are combined together results are shown for uq = 36 GHz. 
Only 39 of 43 pixels for each channel were considered in this analysis. 

^Correlation slope a in units of /xK(MJy/sr)~^. 
'^Temperature anisotropy AT in units of yitK. 

^In this case the best fit n was greater than 15, outside the range of the n grid used. 
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Fig. 1. — The figure shows the OVRO data (open circles connected by sohd lines), the 100 ^m- 
correlated emission (dashed Une), the 12 /xm-correlated emission (dot-dashed hne) and the Ha- 
correlated emission (dotted line), all at 14.5 GHz. 
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Fig. 2. — The figure shows the SP94 Q band coadded data (open circles connected by sohd hnes), 
the 100 /um-correlated emission (dashed hne), the 12 /xm-correlated emission (dot-dashed hne) and 
the i^a-correlated emission (dotted Hne) in the Q band. The vertical line is at the 39th data point 
marking out the 4 data points that have not been used in this analysis (see text). 
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Fig. 3. — (a) 100 /im-correlatcd emission in the Tenerife (open stars, de Oliveira-Costa et al. 2001; 
filled stars, Mukhcrjcc et al. 2001), OVRO (filled triangles), 19 GHz (open triangles), SP94 (open 
circles), Saskatoon (open square) and DMR (filled circles) microwave data, (b) 12 //m-correlated 
emission reported so far in microwave data (symbols are the same as above). The solid curve is 
representative of a spinning dust spectrum, the dashed and the dot-dashed lines represent free-free 
and vibrational dust emission spectrums respectively. 



